Lyddane- Sachs- Teller relationship in linear magnetoelectrics 
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In a linear magnetoelectric the lattice is coupled to electric and magnetic fields: both affect the 
longitudinal-transverse splitting of zone-center optical phonons on equal footing. A response matrix 
relates the macroscopic fields (D, B) to [E, H) at infrared frequencies. It is shown that the response 
matrices at frequencies and oo fulfill a generalized Lyddane-Sachs- Teller relationship. The rhs 
member of such relationship is expressed in terms of weighted averages over the longitudinal and 
transverse excitations of the medium, and assumes a simple form for an harmonic crystal. 

PACS numbers: 75.85.-|-t, 77.22.-d 
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The original Lyddane- Sachs- Teller (LST) relation- 
ship [1| applies to the simple case of a cubic binary crystal 
in the harmonic regime. It relates four macroscopically 
measurable constants as 



e(cxD) 






(1) 



Here e(0) is the static dielectric constant, which includes 
the lattice contribution, e(oo) is the so called "static high 
frequency" (a.k.a "clamped- ion" ) dielectric constant, 
which accounts for the electronic response only, and wl 
(wt) is the zone-center longitudinal (transverse) optical 
frequency [2|. It is remarkable that all microscopic 
parameters (force constants, masses. Born effective 
charges, cell volume) disappear from Eq. ([T]). 

Magnetoelectrics (MEs) are insulators where electric 
fields control magnetization, and conversely magnetic 
fields control polarization; they attracted considerable 
theoretical and technological interest in recent times [3|- 
ll|. The simplest and most studied single-crystal linear 
ME is antiferromagnetic Cr203 [3|, m\m- I^- any linear 
ME the role of the dielectric function e(a;) is played 
by the 2x2 response matrix — called 7^(a;) here — which 
yields the macroscopic fields {D,B) in terms of {E,H) 
at frequency w: 



7^(a;) 






, (2) 



where e, fi, and a are permittivity, magnetic permeabil- 
ity, and ME coupling, respectively. In this Letter we are 
going to show that a generalized LST relationship holds 
in the form 



tr {7^-l(cx))7^(0)}-2 _ loJ 



(3) 



In the simple case of a magnetically inert material (i.e. 
fi = constant, a = 0) the Ihs of Eq. ([3]) equals indeed 
e(0)/£(oo). While in ordinary dielectrics the LO-TO 
splitting is due to the coupling of the ionic displacements 



to macroscopic electric fields, in linear MEs it is due to 
the coupling of both (electric and magnetic) fields on the 
same footing: this is perspicuous in the Ihs of Eq. ([3]). 

The simple form of Eqs. ([1]) and ([3]) requires a 
crystalline system with only a single IR-active mode at 
the zone center. The additional requirement of cubic 
symmetry can be relaxed; it is nonetheless convenient to 
consider only crystals whose symmetry is orthorombic or 
higher; then all crystalline tensors can be simultaneously 
diagonalized (as e.g. in Cr203). This allows us to adopt 
simple scalar-like notations, as in Eqs. ([T]) and ^ with 
the proviso that we deal separately with each principal 
direction. 

Over the years the LST relationship has been extended 
in several ways, to cover cases where more than one IR- 
active mode per direction exists 12|, |l3[, the crystal is 
in a low-symmetry class IJ, |l5| , or even the material is 
noncrystalline and/or anharmonic 
case can be written as 



H-d^. The general 



£(oo) 



(c^2)t ■ 



(4) 



The quantities in the rhs are weighted averages, obtained 
from moments of the appropriate spectral functions. The 
derivation is based on general principles of statistical 
mechanics and does not require an Hamiltonian, even 
less an harmonic one 16|, Il8l |: in the special case of a 
single IR-active harmonic mode per direction Eq. ^ is 
equivalent to Eq. (TTj) . In this Letter we generalize the 
viewpoint of Ref. [18[ to the ME case, showing that 



tr {7^-l(oo)7^(0)}-2 {uj^)i 



2-tr {7^-l(0)7^(oo)} {uj^)^ 



(5) 



where the rhs is defined below, Eqs. (|10l) and (ITT|) . 

The presentation proceeds as follows. We start at a 
very general level without any microscopic assumption 
about the ME medium, and using only very general 
principles in order to arrive at Eq. ([S]). We will then 
apply the general results to a crystalline system in the 



harmonic regime, and finally we will show that Eq. ([S]) 
reduces to Eq. ^ in the single-mode case. 

We write explicitly the linear response matrix of the 
ME medium as the sum of its real and imaginary 
part: TZ{uj) — TZ'{uj) + iTZ"{uj), and analogously for 
its inverse; both TZ{uj) and TZ^^{uj) obey the Kramers- 
Kronig relationships in the form 

7^'H-7^(oo) ^ ^-Tdu'^^^^p^, (6) 



TT Jq w'2 — UJ^ 



From these, it follows immediately that the numerator 
and denominator in the Ihs of Eqs. (j31) and ([5]) are 

tr{7^-l(oo)7^(0)} - 2 = - / — tr{7^-l(oo)7^"H}, 

(8) 

2-tr{7^-l(0)7^(oo)} = / — tr{7^-l"(w)7^(oo)}. 

(9) 
The integrands in the rhs of Eqs. ([5|) and ^ are 
interpreted here as the transverse and longitudinal 
spectral weights, respectively, by means of which we 
define the second moments 



i^'h 



(-^)l 



/o"^c.^tr{7^-lM7^-H} 
/-i^tr{7^-lM7^"H} 

Jo" ^uHv {7^-l'^M7^(oo)} 
/-i^tr{7^-l"H7^(oo)} ■ 



(10) 



(11) 



The reason for the semantics (transverse and longitudi- 
nal) may appear obscure at this point; it will become 
clear when specializing Eqs. ([TU)) and (fTTj) to an harmonic 
crystal — see also Eq. (|26t below. 

In order to arrive at our main result, Eq. ([5]), we exploit 
the "superconvergence" theorem [20[. For large w (i.e. 
for w much larger than all the resonances of the medium) 
TZ"{uj) vanishes; Eqs. ^ and ([7]) yield, to leading order 
in 1/^2, 



■R{Qo)-^n{uj) 



niLu)-^n{oo) 



duj' uj'n{oo)-'-n"iuj') 



(12) 



duj' w'7^-l"(w')7^(oo), 



(13) 

where I is the 2x2 identity. Inversion of Eq. ([T^ to the 
same order gives the alternative expression 

2 f°° 

n{uj)-^n{oo) ^ I + — - / duj' Lu'n{oo)-^n"{uj'). 

TTW^ Jo 



Next we take the trace of Eqs. (|13p and (HH); permuting 
the matrices in the product we get the identity 



duj wti- {7^(oo)"^7^"(w)} 

doj uj tr {TZ{uj)-^"tZ{oo)}. (15) 



In order to arrive at Eq. ([5]) it is enough to put together 

Eqs. (0), (ED, (nni), cn), and i^. 

We stress that at the root of the superconvergence 
identity, Eq. (|15|) . is the assumption that lo = oo actually 
means uj much higher than all the frequencies of ionic 
motions, yet lower than the frequencies of electronic 
excitations [2l- Therefore the clamped- ion response 
7^(oo) is a real matrix. 

Next we address an harmonic crystal. A zone-center 
optical mode is lattice-periodical; it is then expedient 
to consider the energy per cell of a (macroscopically 
homogeneous) solid with a frozen-in phonon distortion. 
This energy is well defined only when a prescription for 
taking the thermodynamic limit is given. We cut a 
sample in the shape of a slab parallel to the principal 
axes; we remind our assumption that all crystal tensors 
are diagonal on them. If the slab is free-standing in 
vacuo, all fields vanish outside {E—D~H—B—0), while 
the value of the fields inside depend on the polarization 
of the mode. Simple electrostatics and magnetostatics 
imply that if the phonon polarization is parallel to the 
slab ("transverse"), both E and H vanish, while if it 
is perpendicular ("longitudinal") D and B vanish t22|. 
The order of the limits (first a slab, then its thickness 
going to infinity) is essential, and the two energies 
(longitudinal and transverse) are indeed different in the 
thermodynamic limit. Similar reasonings apply if the 
lattice-periodical mode is regarded as the k — > limit of 
a finite-k optical phonon [2]. 

If the crystal has N IR-active modes in each principal 
direction, we denote with u;„ the zone-center TO 
frequencies (i.e those with £'=0 and H=0); equivalently, 
o;^ are the eigenvalues of the analytical part of the 
dynamical matrix at k = 0. The free energy per cell in 
function of the normal coordinates w„, E, and H (taken 
as independent variables) is expanded to second order 
as dH 



F{E,H,{u^})=Fo + IY1 



2 2 



- —[e{oo)E^ + 2a{oo)EH + fM{oQ)H^ 

OTT 



(16) 



(14) 



where we adopt atomic Gaussian units 2l| , and il is the 
cell volume. Z* are the mode-effective charges and C* 
their magnetic analogues; notice that in ordinary units 



the normal-mode coordinates would include a factor with 
the dimensions of (mass)^/^, while Z* and (* would 
include a factor with the dimensions of (mass)~^/^. 

The derivatives of F provide the equations of motion 
in the form 



D = - 



AirdF 
n'dE 



= e{oo)E + a{oo)H 



47r Y^ 






in dF 47r ,-^ 



In rx ^ 



dF 



= -U!„U„ 



KE + QH. 



(17) 



We then consider forced oscillations at frequency w, i.e. 
/„ = —Lu'^Un- Elimination of the u„'s from Eq. (J17p 
provides the linear ME response, including the lattice 
contribution; we cast it in compact form as 

47r ■ 



7^'H = 7e(oo) + -^ 



n"iuj) = 



"TT^- 



ZnZl 


(18) 


^l~^' 


[^l-^^], 


(19) 



where the ME lattice coupling vectors are 

-^n =1 ,* I , -^„ — [ '^n ^n ) 






(20) 



Eq. (fT8| is the elegant result obtained in 2008 by J. 
Miguez [5| ; the TO frequencies w„ are clearly the N poles 
oin{uj). 

Actual computations performed for the paradigmatic 
material Cr203 and based on Eq. ([T5)) show that the 
magnetoelectric coupling a(0) is significantly enhanced 
by the lattice contribution [5|]. Clearly, a large coupling 
is the key property to be exploited in device applications 

0- . . 

It is now expedient to express Eq. (|17p in terms of D 

and B. 



-UJ^Un 



TZ-\oo) 






(21) 



A'jr 

— Zln^^ (OO) Y^ Zn'Un' 



zin-\^) 



(22) 



As explained above, in the LO modes D=Q and B={) by 
definition [22]. The first line of Eq. ([^ clearly shows 
that the LO eigenmodes are in general different from the 
TO ones; an explicit N x N diagonalization is needed 
in order to find the LO frequencies. We indicate with 
LUn these frequencies, and we also transform the lattice 
coupling vectors to the LO eigenmodes as TZ~^{oo)Zn -^ 
Zn, in order to write Eq. ((22)) as 



- W Un = -i^n'^n + 



Jt 



(23) 



Elimination of the u„'s from Eqs. (|21l) and (l^5|) provides 
TZ^^iuj) in the form 



^ ^ CjI-uj-^ 



n"-\uo) 



47r2 



Ez„Zt^((i2-c.2). (25) 



Using the above results, we arrive at a very transparent 
expression for the rhs member of Eq. ([S|) : 






(26) 



i.e. {uj'^)t is the weighted harmonic mean of the w^'s, 
with weights -Z^7^~^(oo)-Z„, and (w^)l is the weighted 
harmonic mean of the w^'s, with weights Z^TZ{oo)Zn. 
The two sets of weights are in general different, except 
when the transverse and longitudinal eigenmodes happen 
to be the same. However, "superconvergence" , Eq. (fTSt . 
implies the same normalization in any case: 



Ezt7^-l(^)z„ 



Ezt7^(^)z„ 



(27) 



In the special case where a single IR-active mode exists 
{uj'^)t — ^T and (cl'^)l — cdl'- this concludes the proof of 
Eq. ©. 

A well known bound requires the matrix TZ to be 
positive definite [3|. Since this is based on stability 
arguments, it only concerns 7?.(0) and TZ{oo): the former 
is a genuine static property of the real system, while the 
latter can be regarded as a static property in the infinite 
nuclear mass limit. At any other frequency the matrix 
TZ{uj) accounts for the forced oscillations of the system, 
which is clearly out of equilbrium. Therefore TZ{uj) is not 
required, in general, to be positive definite; because of 
the same reason, e(w) is not a positive real function in 
ordinary dielectrics [2|. 

All of the above results reduce to previously known 
ones for a magnetically inert material, where the 2x2 
matrix TZ{uj) has the unique nontrivial entry s{io). It is 
worth examining Eq. (P5| for Q — ^- 






1 En^jz:)' 



(28) 



The weights in the harmonic means (Z^)^ and (Z*)^ 
are the (squared) transverse and longitudinal effective 
charges, respectively. If (and only if) the transverse 
and longitudinal eigenmodes are the same, then Z* = 
Z*J e{oo). Actually, this is the well known relationship 
between the transverse (a.k.a. Born) end longitudinal 
(a.k.a. Callen) effective charges. 

For an ordinary dielectric the TV poles of e{Ld) are the 
TO frequencies a;„, while the N poles of l/e{uj) are the 



LO ones w„: see Eqs. ([T5)) and ([M]) . Since the response is 
a single-component function, the poles of l/eiuj) coincide 
with the zeros of e(w). Simple considerations about zeros 
and poles of £{uj) eventually lead to the simple expression 






N _ 



(29) 



first found in 1961 by Kurosawa [12|,ll9|. This result does 
not generalize to the ME case, for the good reason that 
the response is a 2 x 2 matrix: the present formulation is 
based on traces throughout, and the inverse of the trace 
bears no simple relationship to the trace of the inverse, 
at variance with the purely electrical case. 

This Letter addresses the linear relationship between 
the pairs [E, H) and (£>, B) throughout, starting with 
Eq. ([2]) onwards. Other pairings are possible. In 
particular the choice [E, B) and (D, H) looks like a more 
natural one for at least two reasons: (i) the microscopic 
forces are determined by the {E,B) pair [2^, and (ii) 
a relativistic formulation can be elegantly cast in terms 
of two four-dimensional field tensors whose entries are 
{E,B) and {D,H), respectively [a, lZ|. Nonetheless, 
the same choice in the present context would not be a 
convenient one. In fact, as explained above (see also Ref. 
[22|). in a transverse mode E — Q and B ^ 0, and in a 
longitudinal one D = and H ^i). 

Throughout this Letter we have stressed the formal 
equivalence of electric and magnetic fields in their 
coupling to the lattice in MEs. However, the orders 
of magnitude of electric and magnetic phenomena in 
condensed matter are not the same. ME effects are 
notoriously small 13], and the corrections to the LST 
relationship in most cases are expected to be small as 
well. In oxides the dielectric constants — either e(cxD) or 
e(0) — are typicall y in the range 2 to 10, while l/x — 1| is 
of the order lO""*^ [23|; in the most studied linear ME, i.e. 
Cr203, even a is of the order lO""*^ [5|, |6[. An accurate 
evaluation of the Ihs of Eqs. ([3]) and ([5]) in "conventional" 
ME materials would require a measurement of all the 
entries of the response matrices Ti. to the same absolute 
error, which could be problematic. More perspicuous 
effects are expected in nonconventional materials [3|, 



24| 



such as those where the ME effect can be tuned 

In conclusion, this Letter shows that the ratio 
e(0)/e(oo) entering the LST relationship must be 
replaced, for a linear ME, by the Ihs of Eqs. ^ and (O, 
whose ingredients are the full ME responses at frequency 
and oo, i.e. static and clamped-ion. In the most general 
case the rhs member of our generalized LST relationship 
is the ratio between spectral moments of the longitudinal 
and transverse excitations of the medium. It assumes 
a simple form for a crystalline ME in the harmonic 



approximation, and finally is identical to the original LST 
one when only one mode is IR active. The relationship 
shows that the LO-TO splitting in a ME originates from 
the coupling of ionic displacements to both electric and 
magnetic macroscopic fields. 
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